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I. MAGAF SYSTEM ADDITIONS 
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SUBROUTINES SUNPK AND $REPK 


SUBROUTINES S$UNPK AND $REPK HAVE BEEN WRITTEN 
RESPECTIVELY TO UNPACK A DATA FRAME FROM RECEIVED 
DATA ORDER TO INSTRUMENT ORDER (DEFINED IN A PREVIOUS 
REPORT) AND TO REPACK THE DATA INTO RECEIVED DATA 
ORDER. SUBROUTINE REPACK IS USED PRIMARILY TO OBTAIN 
DUMPS FOR HARDWARE DIAGNOSTIC PURPOSES. 


A PERMUTATION TABLE (SCYCST THRU #S$CYCND) 
WHICH DEFINES THE PARTICULAR PERMUTATION FOR THIS 
PURPOSE HAS BEEN ADDED TO THE IN CORE SYSTEM, 


USAGE: DIMENSION 1A(245), 18(245) 
CALL SUNPK(1IA,1B) 
CALL SREPK(1A,1B) 
RESTRICTION: JA AND 1B MUST BE DISJOINT OR 
IDENTICAL. 


BECAUSE THE METHOD USED IS GENERALLY APPLICABLE 
TO ANY PERMUTATION AND ITS INVERSE, A BRIEF DIS- 
CUSSION OF METHOD AND A FLOW CHART OF THE PROCEDURE 
WILL BE GIVEN. 


THEORY: 

ANY PERMUTATION CAN BE REPRESENTED AS A PRODUCT 
OF DISJOINT CYCLES, WHERE A CYCLE IS MERELY A TRACING 
OF ELEMENTS WHICH PERMUTE INTO EACH OTHER. (SEE ANY 
TEXT ON MODERN ALGEBRA). FOR EXAMPLE, USING PERMUTA- 
TIONS OF A SQUARE, ROTATION THROUGH 90°: 


4 | 


CAN BE REPRESENTED AS (1,2,3,4) i.e. | GOES TO 2 
GOES TO 3 GOES TO 4 GOES TO 1 AND REFLECTION ABOUT 
A DIAGONAL 

hin 


4 3 2 
CAN BE REPRESENTED AS (1) (3) (2,4). 


THE ADVANTAGE OF THIS REPRESENTATION IS THAT THE 
INVERSE PERMUTATION IS IMPLICITLY DEFINED BY THE PER- 
MUTATION ITSELF - READ THE CYCLES BACKWARDS. 


METHOD: 

AN INDEX TABLE 1S PREPARED WITH NEGATIVE TERMI- 
NATORS AT BOTH ENDS CONTAINING THE CYCLES AUGMENTED 
BY THE START INDEX OF EACH CYCLE. 

FOR EXAMPLE THE TABLES FOR THE ROTATION AND 
REFLECTION WOULD BE: 

ROTATION REFLECTION 
=] -1 
| 


2 
3 
4 
! 
| 


THE PROGRAMS THEN SCAN THE CYCLES IN THE INDEX 
TABLE IN THE PROPER DIRECTION TO PERFORM THE PER- 
MUTATION. 


FLOW CHART: 


C oe) 


GET FIRST 
(LAST) INDEX 
ENTRY 


GET NEXT 
INDEX ENTRY 


<foo> (RETURN 
ee 


STORE FIRST 
CYCLE INDEX 


GET CURRENT 
VARIABLE 


STORE 
VARIABLE 
{INTO PREVIOUS 
ACATED SPACE 


GET NEXT 
INDEX 


~~ NEXT 
INDEX = FIRST 
YCLE INDEX 


NO 


STORE CURRENT 
VARIABLE 
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NAME SCYCST, BCYCND 
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DATA STORAGE AND RETRIEVAL SUBROUTINES 


Tne DATA STORAGE AND RETRIEVAL SUBROUTINES ARE DESIGNED TO 
FACILITATE THE STORAGE AND RETRIEVAL OF TIME SERIES DATA 
WITHOUT REQUIRING INORDINATE AMOUNTS OF CORE. THERE ARE 
S'X SUBROUTINES, EACH OF WHICH WILL BE DESCRIBED. 


ROUTINE: OPSN 
CALLING SEQUENCE: CALL OPSN (BUFFER,NVAR) 


BUFFER - FLOATING POINT ARRAY OF DIMENSION AT 
LEAST 60*NVAR (120*NVAR FOR FIXED POINT) 


NVAR - NUMBER OF VARIABLES 


OPSN - MUST BE CALLED BEFORE ANY OTHER ROUTINES TO 
OPEN THE DISKFILE (PLOTFL ON UNTIL 30) AND 
SETUP FOR OTHER ROUTINES 


ROUTINE: | ADDCMP (ADD COMPONENT) 
CALLING SEQUENCE: CALL ADODCMP(A,1) 
A - VARIABLE TO BE ADDED 
| - INDEX OF VARIABLE SNNAR 


ROUTINE: RETRIV 
CALLING SEQUENCE: CALL RETRIV(VECT,1#) 
VECT - VECTOR WHICH IS TO BE RETURNED 


| - INDEX OF VARIABLE 


ROUTINE: CLOSE 
CALLING SEQUENCE: CALL CLOSE 


CLOSES DISKFILE AND UPDATES IT. REOPENS IT 
FOR LATER USE 


NOTE: FOR PROPER PROGRAM FUNCTIONING, ADDCMP MUST 
HAVE BEEN CALLED THE SAME NUMBER OF TIMES 
FOR EACH INDEX, AND THAT NUMBER MUST BE A 
MULTIPLE OF 60. 


tat. sateen ts lla 


RR RE oe we 
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AOUTINE: CLOSE 
CALLING SEQUENCE: CALL CLOSE 


CLOSE DISK FILE AND UP DATE IT, REOPEN IT 
FOR LATER USE. 


NOTE: FOR PROPER PROGRAM FUNCTIONING, ADDCMP MUST HAVE 
BEEN CALLED THE SAME NUMBER OF TIMES FOR EACH 
HNDEX AND THAT NUMBER MUST BE A MULTIPLE OF 60. 


ROUTINE: OPSNI, CLO$) 


CALLING SEQUENCE: CALL (OPSN!) (TFIRST, ILAST) 
(CLO$s!) 


(FROM) pis, 


(OPSNI RETRIEVES) garg getween LFIRST AND (LAST i } 


(CLO$T STORES ) 


EXAMPLE: A TAPE EXISTS IN THE FOLLOWING CARD IMAGE FORMAT. 
RECORD | - Alphabetic (0. 


RECORD 2 - ny - Xp - XzQ - OBSERVATIONS OF 20 
VARIABLES (20F4.9) n<2401 


PROBLEM: TO STORE THE DATA MATRIX ON DISK, LATER 
TO PRINT OUT THE DATA VARIABLE BY VARIABLE. 


JIMENSTON T0040) THUPRE RE LA00) Xe? 
COMMON 4XRRKZ 1LOON 
REWIND @t 


READ ChLetOO) TO 
100 FORMAT (@0A2) 
Nev 


GALL OPEN HUFFE Re 20) 
A READ Cebelold x 
POL FORMAT C20OF 40) 


TPC TOCHR Ge LD eee 
2 UACKSRPALCE Al 
“oO TO 1 


3 Nene} 
DO S Lele @o 
S CALL ADUCMPECXCLDe Tb) 


GO TO )1 
& NADD SEMOD INI6O) 
IF (NADDe EQe 0) GO TO YO 


OO & T® beNnaAdvD 
DO © Jetetd 
© FALL ADLCMP(Oe0eU) 


YO CALL CLOSI Coe O(al)) 
STOP 
END 
UTMENSTON LOCa@tie VECT R400) 
CALL OPENIVECTE@O) 
CALL OP ONT ETD 1D14)) 
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Eau 
EQu 
EQu 
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INCR 
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ADL 
STA 
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STA 
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OP$N, CLOSE ,RETRIV,ADDCMP,OP tr 
#SE Chose} 
422 

1 

2 


&SE,2,4,0 


SAVE 
RDW,B 
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CnLL 
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RETA 


SAVE 
PECHA 
FORt3 
FOE, 36, .i 


RETAN 


CLOSE 


SAVE 
NOOMP 
Noone 
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QQOaZO RETRN 800068 
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OO@1SL WHERE OQdea3e 
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A 
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RETL 
SAVE 
STX 
WORK 
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CORe1LSS 
QOS AT 
OOM 0 
OOOIE 4 
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ATU On rons 
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CLORE 
TCOMP 
OP AN 
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STA 
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OOeB21a 
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Eats 
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VT 
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Cc. 
TEKTRONIX PLOTTING ROUTINES 


A SERIES OF PLOTTING ROUTINES HAS BEEN ADDED TO THE SYSTEM 
TO ENABLE PLOT FILES IN VARIAN DATAPLOT FORMAT TO BE OUTPUT 
TO THE TEKTRONIX 4O14&. THESE ROUTINES SCALE A STATOS PLOT 
SO THAT ONE INCH ON THE STATOS EQUALS ONE INCH ON CRT, NOT 
ON THE HARDCOPY. THESE ROUTINES ARE TRANSPARENT TO THE 
USER, WITH THE FOLLOWING EXCEPTIONS: 


1) CARE SHOULD BE USED WITH NEGATIVE ORIGINS (in inches) 


2) PRINTS BEYOND THE SCREEN LIMITS END UP AT THE RIGHT 
HAND AND TOP LIMITS OF THE SCREEN 


3) ALL CHARACTERS ARE PLOTTED AT THE CURRENT 
CHARACTER SIZE (see TEKFNC below) AND UPRIGHT-ORIENTATION 


4) SPECItAL CHARACTERS HAVE NOT BEEN IMPLEMENTED 


5) MINIMIZATION OF STATOS SORT AND PLOT TIMES MAY 
INCREASE TEKTRONIX PLOT TIMES 


IN ADDITION, A SUBROUTINE TEKFNC HAS BEEN WRITTEN TO ENABLE 
THE USER ACCESS TO TEKTRONIX FUNCTIONS. 


USAGE: 


TO OUTPUT TO THE TEKTRONIX, ADD THE FOLLOWING SUBROUTINES 
TO ANY DATAPLOT PROGRAM: 


SUBROUTINE DPSORT 
CALL CRTIPLT 
RETURN 

END 

SUBROUTINE DPPLOT 
RETURN 

END 


TO USE TEKTRONIX FUNCTIONS: 


CALL TEKFNC(1) 


FUNCTION 
LARGEST CHARACTERS 
2ND LARGEST CHARACTERS 
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(SMALLEST) CHARACTERS 
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QOAQOeL 
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QOORBO4A 
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088006 
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eedae? 
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OQOO4L 
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BOON44 
OBOO4S 
QOOOAG 
BOdO047 
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OOHOE1 
eoedee 
000423 
QAQOOD 


Oa6S05 
QOOOWAA 
1O0000 
QOL410 
QAQ133 
QeAHadYD 
QOOOOO 


QooOGO1L? 
QIOHO OHO 
OOLO1LE 
POOLS 


eodBSOS 
Haqvaes 
1OvVOOS 
O1800e02 
ONO LAS 
ooo 
OOLKBOO 


OOBSOS 
OoOvVOLS 
LQModoOo 
OoOodLO 
@OOLZS 
QOOOIOO 
QAQAAD 
OL4tat 
O54005 
OOGOAA 
@00836 
eacgzee 


OVAVaa 
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HOOOOA 


OOHAL?T 
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CRTPLT ENTR 
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READ 
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DBR DER 
EXT 
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Jn 
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OQOaSE 
ANOS? 
YOO 
YOOOH IL 
AAA 
QMO S 
AOI 4 
OOS 
OMA 
OVOAG? 
YAIOQTVA 
YOM? 
VYYaaTvea 
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adOLed 
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@BOOLEL 
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OdoOLea+4 
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OoON1LA? 
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YOALAL 


OdaLAe 
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L204e3 
QYGASG? 
OODOAL 
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OVALE 
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Oey every 
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OO ery 
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DATA ANALYSIS 


I. 


A program to produce magnetograms from the AFGL 
Magnetometer Network data was written. The program is 
to be run on the Varian computer and the magnetograms 
are hard copies of a Cathod Ray Tube (CRT) display. 


The program is run by typing in on the CRT keyboard 
/LOAD, DEMO 

The program DEMO then prompts the programmer 
ENTER START HOUR OR 99 FOR BOT. 
93 FOR CURRENT TAPE POSITION, 97 FOR DISK DATA 
If you want to begin the magnetogram at a particular hour 
enter that hour as hh. If the magnetogram is to begin 
at the beginning of the TAPE (BOT) enter 99, if it is to 
begin from the current position of the tape enter 93. 
If the magnetogram is to be made from data already on 
disk enter 97. 
ENTER START TIME FOR GRAPH IF DIFFERENT 
The magnetograms should begin on the hour, but if the data 
to be plotted do not begin on the hour, enter a starting 
hour. For example say the data begin at 17:37, the 
programmer may enter 17 as the starting hour of the 
magnetogram. otherwise enter RETURN 
ENTER START DAY IF DIFFERENT FROM CURRENT TAPE DAY 
If the tape is positioned on day 213 for example and you want 
a magnetogram beginning on day 214 then enter 214. Otherwise 
enter RETURN. 
ENTER NUMBER OF HOURS AND NUMBER OF PERIODS IF NOT 1 (HHPP) 
Enter the length of the magnetograms in hours and the 
desired number of magnetograms. If you want seven eight 
hour magnetograms enter 0708. If only eight hour magnetogram 
is wanted enter 08. 
AUTO HARD COPY (Y,N)? 
If you want hard copies of the magnetograms produced auto- 
matically enter Y, if not enter N. 
INDIVIDUAL PLOTS, OVERLAY PLOT, OR BOTH (I,0,8)? 
If you want an individual plot of each station and no overlay 
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enter I, if you want an overlay and no individual plots 
enter 0. If you want both enter B. 


DEMO is a highly automated program. Once all the 
entries have been made all that is required of the 
programmer is that he change the tapes if the magnetogram 
requires more than one tape and that he terminate the 
program. The program is terminated by simply entering a 
/ RETURN. 


An example of the prompting sequence and the resulting 
magnetograms is shown in Figures 3 thru 9. These are the 
plots of X, Y, and Z components of the seven network 
stations, in crder: MA, FL, MI, WI, SD, CA, WA. Figure 10 
is a composite of all seven stations. 


Figures 1 and 2 are plots of LINCHK for the selected 
date of the magnetograms, and give an overall view (for 
that date) of the quality of operation of the network for 
each station. 


To use LINCHK 
Load the tape of interest and type on the CRT 
keyboard /LOAD,LINCHK 
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